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Within the framework of a Modified Poisson-Boltzmann theory accounting for steric effects of 
microions, we prove analytically that the effective pair interactions between like-charge colloids 
immersed in a confined electrolyte are repulsive. Our approach encompasses and extends previously 
• known results to the case of complete confinement, and further incorporates the finite size of the 

microions which is absent in the standard Poisson-Boltzmann theory. 
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! PACS numbers: 05.70.Np, 64.10+h, 82.60Lf, 82.70.Dd 

Recent experiments show convincingly the existence of long-range electrostatic attractions between like-charge 
colloids immersed in an electrolyte, in particular in the vicinity of a charged wall or when the particles are confined 
in a slit These striking observations are inconsistent with the well established theory of Derjaguin, Landau, 

. Verwey and Overbeek (DLVO) fl and constitute an important controversy in colloid science. Inasmuch as direct 
measurements of the effective pair interactions in the bulk of a suspension confirm the validity of DLVO theory with 
repulsive interactions confinement effects are thought to play a major role and have been included in numerical 
H and theoretical |6) investigations of the electrostatic forces between colloidal spheres. These studies reported a 
range of inter-particle distances for which the computed forces were attractive. However, Neu & rigorously proved 
that the Poisson-Boltzmann (PB) model considered in || necessarily leads to repulsion. His proof applies within PB 
theory (which forms the basis of the DLVO potential) in the specific case of Dirichlet boundary conditions (constant 
electrostatic potential ip on the confining surface, the experimental relevance of which is dubious). Sader and Chan 
I ' extended the PB argument to a broader class of boundary conditions || . 
, The PB theory is a continuum mean-field approach neglecting statistical correlations between microions, that are 
assumed to be point-like. It is known to grossly overestimate counter-ion concentrations close to a charged surface (eg 
CN| ! that of a polyion). A Modified Poisson-Boltzmann (MPB) theory including steric effects was put forward to overcome 
this shortcoming This approach retains the simplicity of the original PB theory and is more tractable than 

the other attempts to improve upon PB by inclusion of steric repulsion p3| , p^I[ . 

In this article, we consider the MPB model in the situation of a mixture of microions (with overall global elec- 
troneutrality). Within this framework, we analyze the interactions between a pair of like-charge polyions immersed 
in a confined electrolyte of permittivity s. The colloids may be of arbitrary shape and the confining region TZ is a 
cylinder of arbitrary cross-section. As in |7|||, the only requirement is that the electrostatic potential ip possesses 
mirror symmetry with respect to the mid plane z = between the colloids. Unlike in |7j|^] where TZ is of infinite lateral 
extension, we allow the confining region to be of finite volume (situation hereafter refcred as finite confinement, which 
can be that of a pore or of a closed Wigner-Seitz cell). The geometry considered here encompasses many cases of 
H , experimental relevance, such as two colloids confined in a slit, in a charged pore or in the vicinity of a wall We 
consider boundary conditions where the medium TZ' outside TZ is a dielectric continuum of permittivity r with the 
possibility of a uniform density of surface charges a on the boundary dTZ (as in ^|). This class of boundary conditions 
is large, as for a = the limit e' — » reduces to Neuman boundary conditions with a vanishing normal electric field on 
dTZ, whereas e' — > oo corresponds to the standard Dirichlet condition on dTZ. We show rigorously that the inclusion 
of excluded volume effects within the MPB approach does not change the sign of the force which remains repulsive 
as in PB. As the MPB gives back the PB theory in the well behaved limit where steric effects disappear, our results 
also include those reported in J^,^|, and extend them to the situation of finite confinement. 

The MPB theory can be obtained from the free energy functional T — Tq ou \ + T cn t of the inhomogeneous fluid of 
microions contained in TZ. Given N species a with a local density c a (r) and charge number z a , the Coulombic energy 
contribution to T reads 
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T Co ^ = \ f Pc {r)G(r,r')p c {r')drdr>, (I) 
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where p c (r) — ^ Q z Q ec a (r) is the local charge density and G(r,r') is the Green's function inverting the Laplacian 
inside TZ with the required boundary conditions. In the absence of confinement, G reduces to the Coulomb potential 
l/[47re|r — r'\]. The entropic contribution is 
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where A Q are irrelevant length scales and kT = /3 _1 is the thermal energy. For the sake of simplicity, all types of 
microions are assumed to have the same size a. The last term in Eq. (|^) mimics the entropy of the solvent molecules 
while the remaining ones are the ideal entropies of the mixture of co- and counter-ions. In the limit where a — > 0, the 
steric correction disappears and the classical PB expression is recovered. The equilibrium density profiles are those 
which minimize the free energy functional T subject to the constraint that 



c a (r)dr = N a 



■r 



with N a the total number of ions a inside 1Z. The minimization yields: 

Pa exp(-/3ez a V) 
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where ip(r) = j n p c (r')G(r,r') dr' . Since the potential ip can be shifted by an arbitrary constant, the prefactors 
p a have individually no physical significance. In the limit where 1Z is of infinite extension, ip is conveniently chosen 



to vanish far from the polyions and we get p a 



„bulk 
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In the case of complete confinement, two 



experimental situations can be distinguished: a) the solution is of fixed ionic composition (canonical description) and 
the prefactors p a are determined by the normalization constraint ([5j) , or b) the dispersion is in osmotic equilibrium 
with a salt reservoir, treated within MPB for consistency. We then have p a — c a /(l — ^ Q a 3 c a ), where c r a denotes the 
concentration of species a in the reservoir (subject to the electroneutrality constraint J2 a z a c a — 0). As expected, 
the density profiles tend to their PB counterparts as the excluded volume a 3 —> 0. 

The force acting on a colloidal particle follows from integration of the stress tensor II over the surface S of the 
polyion (see Fig 1): 



F = j> U-ndS, 
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where n is the unit vector pointing outwards from the surface of integration. The stress tensor can be obtained by 
considering the mechanical equilibrium condition of a fluid element of microions: the balance between the electric 
force and the osmotic constraint can be written 



-VP + p c E = 

<^ d^P = -p c , 



from which we deduce the expression of the osmotic pressure: 

p (tp) = ~r ln 1 + / ] fl3 P" exp(-/3ez Q V) 
Rewriting Eq. (0) as V ■ II = finally yields 
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where E = —Wijj is the electrostatic field, D = eE and J denotes the identity tensor. It is worthwhile to mention 
that Poisson's equation can be written in the form 
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As the mechanical equilibrium condition invoked above implies that II is divergence free, the surface of integration 
in can be deformed to any surface S' enclosing colloid S only, as stressed in It is convenient to choose 

S' = S U S U El (see Fig. 1). The z-component of the force F then reads 
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On the cross sections So and Sj,, z ■ II • w = [-Dx-E^ — D ■ E/2 — P(ip)]z ■ n, while on the lateral surface £, 
z ■ II ■ n = —<tE z + D' n E' z . In the previous expression, primed symbols refer to the region TV outside TZ (with 
permittivity e') and the dielectric boundary conditions have been used (D' n — D n = a and E z = E' z ). The last term 
in Eq. (hlh can be recast by considering the divergence free tensor II' in TV obeying similar constitutive relations as 
(^) and (jg) (we allow the continuum in TV to contain an electrolyte solution; on the other hand, if no salt is present 
outside TZ, the corresponding pressure P(ip') vanishes in (g)). We get 



D' n E' z dS = 
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where S U £ is the Oxy plane. 

Upon substitution of (Oh into (|TT|) , the axial force can be expressed in the form 
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where primes have been omitted for the part of Oxy belonging to TV (ie S ) : when unambiguous, the same notations 
are hereafter used for the fields inside and outside TZ, with e standing for the permittivity e in TZ and for s' in TV . 
In the last term of (|l3|), (E z ) z =o = due to the mirror symmetry and (E z ) z =l vanishes only in the limit L — > oo 
considered in |^,|| . Expression ( |l3| ) is conveniently recast invoking the identity 



[(£> ■ E) z=0 - (D ■ E) Z=L ] = e (E z=0 - E z=L f + 2 E Z=L ■ {D z=0 - D Z=L ) 



and the relation 
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Equation. (O) follows from a standard Green identity. The line integral in (H4) has been reexpressed remembering 



the dielectric boundary condition D' n — D n = a and that J as (ip z =L — ^pz=o) dl — J s E z dS. Expression ( |T5[ ) was then 
obtained making use of Poisson's equation (|l0|). 
Gathering results, we finally obtain: 
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It can be checked that since the coefficients p a are positive, the osmotic pressure P is a convex-up function of its 
argument ip. Consequently, the first integral in ( jlif ) is positive from which we conclude that F z > 0. This rigorous 
result holds irrespective of the specific boundary conditions to be applied on the colloids and is independent of the 
sign of the surface charge a. The complete confinement of the electrolyte solution in TZ results in an enhanced 
repulsion with respect to the geometries considered in for which the last term of ([l6]) vanishes. As stressed by 
Neu 0, the non convexity of the pressure P(ip) is a necessary condition for attractive interactions. Note that this 
result holds beyond PB or MPB theories, in fact in any continuum description of the electrolyte solution relying 
on the local density approximation: once such a theory provides the functional dependence of the charge density p c 
on the electrostatic potential ip (or equivalently of the pressure P appearing in the stress tensor (^)), d^,p c < (or 
equivalently d\P > 0) is a sufficient condition for repulsive pair interactions. 
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While the numerical results reported in |3| appear erroneous, the reason for the attraction found in M could lie in the 
effective charge assigned to the polyions to account for the presence of tightly bound counterions and polyion-microion 
excluded volume. Our alternative approach to incorporate steric effects gives rise to repulsion. A more definite answer 
could be obtained by considering the generic density functional of Biben et al. Jl5| for electric double layers which 
allows to cope with the finite size of the microions (treated as charged hard spheres) and the molecular nature of the 
solvent (considered to be a mixture of dipolar hard spheres). Our results nevertheless show the robustness of repulsive 
interactions within mean-field theories and may point to the importance of correlated microion density fluctuations 
(neglected by density functional theories) in interpreting the experimental data. 
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Figure 1 (Trizac and Raimbault) 




FIG. 1. Schematic side view of the geometry considered. The closed cylinder S' is made of the lateral surface E and the 
two cross sections perpendicular to the Oz axis (Eo and El). The thick line denotes the boundary dlZ. 
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